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We study asymptotic properties of diffusion and other transport processes (including self-avoiding 
walks and electrical conduction) on large randomly branched polymers using renormalized dynamical 
field theory. We focus on the swollen phase and the collapse transition, where loops in the polymers 
are irrelevant. Here the asymptotic statistics of the polymers is that of lattice trees, and diffusion 
on them is reminiscent of the climbing of a monkey on a tree. We calculate a set of universal scaling 
exponents including the diffusion exponent and the fractal dimension of the minimal path to 2-loop 
order and, where available, compare them to numerical results. 



PACS numbers: 64.60.ae, 05.40.-a, 64.60.Ht 



I. INTRODUCTION 

In 1976, de Gennes coined the nickname ant in the 
labyrinth for the problem of random walks on a fractal 
structure such as a percolation cluster near the percola- 
tion point [jjL In 1982, Wilkc, Gefen, Ilkovic, Aharony, 
and Stauffer 2| introduced the parasite problem as a vari- 
ant of the former problem in which a random walk takes 
place on a lattice animal, i.e., a large percolation clus- 
ter in the region right below the percolation point. The 
statistics of lattice animals belong to the same universal- 
ity classes as the statistics of randomly branched poly- 
mers (RBPs) in dilute solutions in their swollen phase 
and at the collapse transition (9-line), respectively 0- 
[f| . Because cycles or loops of the animals and RBPs are 
irrelevant for these universality classes, they are the same 
as for lattice trees. In other words, the animal or plant 
on which a random walk takes place in these universal- 
ity classes is tree-like. Thus, the random walks we are 
studying in this paper remind us more of a monkey on 
a tree than a parasite on a loop-containing animal, and 
we prefer to glean our nickname from primatology rather 
than parasitology. 

Topologically, trees are 1-dimcnsional structures. 
Hence, the backbone between two separated points on a 
tree consists of an single unique self avoiding walk (SAW). 
Therefore, all the scaling dimensions d a of the backbone 
fractals - the backbone itself, the minimal, the mean, 
and the maximal path, the electrical resistance, the red 
bonds, etc. - are identical: 



d.B = rfmin = ^SAW = ^max = dli — d 



red ■ 



(1.1) 



A clever monkey will climb the minimal path with frac- 
tal dimension d m i n to get a coconut at the end of the 
path. This is in contrast to diffusion, which corresponds 
to an erratic motion of a monkey. Diffusion, on a frac- 
tal medium with fractal dimension df is described by the 



scaling law 



((v(t)-v(0)Y) = R%f(t/R d N -). 



(1.2) 



where r(£) is the position of the walker (here, the monkey) 
at time t, Rn ~ N 1 ^* is the gyration radius of the fractal 
with mass (number of sites) N and / is a scaling function 
with the properties 



/(*) 



1 for x — > oo 

x 2 / dm for x — >• 



(1.3) 



As a result of Einstein's relation d w = df + d_R for the 
fractal dimension of the random walk, it follows that 0,0] 



dw = d f + d n 



(1.4) 



As mentioned above, the monkey on a tree problem 
has been around under a different name for some 30 
years now. For background, we refer to the review arti- 
cle on diffusion in disordered media by Havlin and Bcn- 
Avraham p}. In recent years, significant advancements 
have been made in numerical simulations on problems 
different from but closely related to the monkey on a 
tree problem. A sophisticated Monte Carlo algorithm 
has been used to simulate lattice animals and trees in 2 
to 9 dimensions [8[ and to measure their static scaling 
exponents with high precision [9j,[l0(. Furthermore, sim- 
ulations have been performed to determine with high pre- 
cision the fractal and multifractal dimensions of SAWs on 
percolations clusters in 2 to 4 dimensions [ll|,[l2|]- Hence, 
we feel that state of the art simulations of diffusion and 
transport and lattice animals in dimensions suitable for 
reliable comparison to field theory have become within 
reach. Thus, we think it is worthwhile to take a fresh 
look at the monkey on a tree problem with field theoretic 
methods. 

The static fractal dimension df = 1/up of the ran- 
domly branched polymers or trees are well known 0, [f| . 
Here, we apply renormalized dynamical field theory 



to calculate d m i n and the related exponents in an e- 
cxpansion to 2-loop order. Since d m j n is equal to the 
dynamical exponent z, the scaling exponent of the time 
a monkey needs to reach a coconut on a tree, we can 
and will calculate d m i n via calculating z of a stochastic 
process that generates RBPs [5|. 



II. MODEL AND FIELD THEORETIC 
APPROACH 

This section serves 2 purposes. First, we review the 
field theoretic model for RBPs that we have developed 
recently [5j . Some of the steps involved in its derivation 
are far from trivial, and its symmetry contents is rich 
and interesting. Hence, we think it is worthwhile to re- 
view the model in some detail. This will also have the 
benefit of making the present paper more self-contained. 
Second, we discuss in broad terms the diagrammatics re- 
sulting from our dynamical model for the swollen phase 
and the collapse transition. We will make the observa- 
tion that the dynamical self-energy diagrams decompose 
into a quasi-static part and a SAW-part that contains 
all the frequency dependence. This observation is a key 
to the subsequent sections as it simplifies the dynami- 
cal field theory for the transport and diffusion exponents 
considerably. 



A. Creation of randomly branched polymers 

Our field theoretic model for RBPs is based on the 
idea of generating their statistics through a mesoscopic 
stochastic growth process. It is well known that the 
general epidemic process (GEP) [13( leads to random 
structures with the properties of percolation clusters [14l — 
Il7| which are, depending on the parameter- values of the 
GEP, below, at, or above the percolation point. The pri- 
mary density-fields describing this process are the field 
of agents n(r,t) and the field of the inactive debris 
m(r,t) = Xj_ dt'n(r,t') which ultimately forms the 
polymer cluster. The following extension of the GEP is a 
modification of a process that we have introduced for the 
description of tricritical isotropic percolation [18| . The 
non-Markovian Langcvin equation describing such this 
process (or rather its universality class) is given by 

pi 
\~ 1 d t n= (\/ 2 + cVm ■ v] n - (r p + g'm + ^—m 2 \n + ( . 

(2.1) 
Here, the parameter r p tunes the "distance" to the perco- 
lation threshold. Below this threshold, in the absorbing 
phase, r p is positive which wc assume throughout this 
paper. In this case the typical final clusters of the de- 
bris generated from a source q5(r)S(t) of agents consists 
of N — (J d d r m(r,oo)) w q/r p particles of the debris, 
and has a mean diameter (gyration radius) ~ \j Jr^. 
However, here we are not interested in these typical clus- 
ters. Rather, we are interested in the large non-typical 



clusters, the rare events of the stochastic process, with 
N 3> q/r p (in this sense q is a small parameter) . We know 
from percolation theory [19| that these clusters belong to 
the universality class of lattice animals. Hence, they are 
the same in a statistical sense as randomly branched poly- 
mers (RBPs) as far as their universal properties go. The 
gradient-term proportional to c describes the attractive 
influence of the debris on the agents if c is negative (as 
a negative contribution to g' does). In principle, other 
gradient-terms like mV 2 n and nV 2 m could be added to 
the Langevin equation. However as long as we have any 
one of these gradient terms into our theory, an omission 
of the others has no effect on the final results, and we 
choose to work with the term proportional to c only for 
simplicity. For usual percolation problems (ordinary or 
tricritical), all of these gradient terms are irrelevant. As 
long as g' > 0, the third order term f'm 2 n is irrelevant 
near the transition point and the process models ordinary 
percolation near r p = [151 ] or non-typical very large clus- 
ters, the swollen RBPs, for r p > 0. We permit both signs 
of g' (negative values of g' correspond to an attraction of 
the agents by the debris, see above). Hence, our model 
allows for a tricritical instability (tricritical percolation 
near r p = fig or the collapse transition of the RBPs 
for r p > @). Consequently we need the third order 
term with /' > (representing self-avoidance) to limit 
the density to finite values in these cases. Physically it 
originates from the suppression of agents by the debris. 
The Gaussian noise-source £(r, i) has correlations 

CM)C(r',£') = (A- V(r, t)6(t - t') - fn(r, t)n(v' , *')) 

Xd(r-r'). (2.2) 



The first part of the noise correlation (|2.2[) takes into 
account that the agents can decay spontaneously, and 
thus g > 0. The term proportional to / simulates the 
anticorrclating or correlating (from attraction) behavior 
of the noise in regions where debris has already been 
produced. If the coupling constant / becomes negative, 
attraction effects prevail. For ordinary percolation this 
term is irrelevant. The form of the dependence of both 
parts of the noise correlation on the field is mandated by 
the fact that the process has to be strongly absorbing in 
order to model RBPs. 

In the past, there have been misconceptions about the 
relation between absorptivity and the form of the noise 
correlation, and we think that it is worthwhile to address 
this relation here in a little more detail. As we just men- 
tioned, our RBP-generating process is like all percolation 
processes strongly absorbing, i.e., its extinction proba- 
bility (probability that the process becomes extinct in a 
finite time-interval (0, t)) is larger than zero, P ex t(i) > 
for t < oo. To guarantee this property, all terms in the 
equation of motion (|2.1[) must contain of course at least 
one power of n. Moreover, and this is the important 
point here, the expansion of the noise-correlation func- 
tion (12. 2p in the density n must begin with a linear term 
[2CJ | . With respect to enforcing strong absorptivity, a 



quadratic part in n of the noise correlation is insufficient. 
The upshot is that phenomcnological considerations are 
sufficient to uniquely determine the form of the noise. 
There is no need to resort to a microscopic formulation 
in terms of master equations or the like to figure out the 
proper relevant contributions to the noise correlations for 
the current process nor is it for percolation processes in 
general. 



B. Field theoretic functionals 

To proceed towards a field theoretic model (for the 
general method of field theory in statistical physics see, 
e -<?-i 1211123), the Langevin equations (|2.fj) and (|2.2[) are 
now transformed into a stochastic response functional in 
the Ito-sense [H HHH 



JeGEp[h,n] = d d x{\ I dth(\- l d t -V 



cVm ■ V + r, p + g'm -\ m n ) n 



— I A / dt hn 



n 



(2.3) 



where n(r, t) denotes the (imaginary) response field con- 
jugated to n(r,t). With this functional, we now have 
a vantage point for the calculation of statistical quan- 
tities via path-integrals with the exponential weight 
cxp(— J7oGEp)- When a source-term (j, ft) is added, where 
j(r, t) = qS(r)S(i) and (..,..) denotes an integral of a 
product of two fields over space and time, this functional 
describes, in particular, the statistics of clusters of debris 
generated by the stochastic growth process (|2.1j) from a 
source of q agents at the point r = at time zero. Denot- 
ing by Tr [. . .] the functional integration over the fields 
with boundary conditions n(r, —00) = n(r, 00) = 0, we 
generally have for the generating functional 

Z\j,j]=Tr[exp(-J eGEF [n,n]+(ln)+(j,n)] =1 (2.4) 

if the arbitrary sources j or j are zero. The first property 
follows from causality whereas the second one originates 
from the absorptive properties of the process. Note that 
the role of causality and adsorptivity can be interchanged 
by the duality transformation m(r, t) < — > —n(r, —t) |15l . 

MM- 

Averaging an observable O [n] over final clusters of de- 
bris (the RBPs) of a given mass N generated from the 
particular source j(r, t) — qS(r)S(t) leads to the quantity 

El nam 



(0) N V(N) = (0[n]6(N - M) cxp((j, n)) ) 



cGEP 



Tr 



0[n]5(N - M) cxp(- Jcgep + qn(0, 0)) 



qTv [O[n]n(0, 0)S(N - M) <sq>(-J eG Ep) 



where 



V(N) = (S(N - M) exp(gn(0, 0))) 



cGEP 



(2-5) 



(2.6) 



is the probability distribution for finding a cluster (a 
RBP) of mass N. 



M= d d rdt\n{r,t) = 



d''i 



>(r) 



(2.7) 



is the total mass of the debris. The field m oc (r) = 
m(r, t = 00) describes the distribution of the debris af- 
ter the growth process terminated. Since the probabil- 
ity distribution should be proportional to the number of 
different configurations, we expect by virtue of universal- 
ity arguments the following proportionality between the 
probability distribution V(N) and the animal number 
An for asymptotically large N: 



V(N) 



Nk q n A n - iV 1 " 



Po 



(2.8) 



where kq and po are non-universal in contrast to the uni- 
versal "entropic" scaling exponent 6. The factor N in 
Eq. (|2.8[) arises because the generated clusters are rooted 
at the source at the point r = 0, and each site of a given 
animal may be the root of given cluster. 

In actual calculations, the delta function appearing 
in averages like in Eqs. (|2.5[) and (|2.6p is hard to han- 
dle. This problem can be simplified by using Laplace- 
transformed observables which are functions of a vari- 
able conjugate to N, say z, and applying inverse Laplace 
transformation (where all the singularities of the inte- 
grand lie to the left of the integration path) in the end. 
The switch to Laplace-transformed observables can be 
done in a pragmatic way by augmenting the original 
JeGEP with a term zM. and then working with the new 
response functional 



Jz = JcGEP + ZM 

As an example, let us consider 



(2.9) 



(0) N V(N)= |ie^ 

x (O[n}cxp(-zM+qh(Q,0))) cGEp . (2.10) 

Note that the relationship between V{N) and An given 
in Eq. (|2.8|) signals the existence of a singularity ~ (z — 
z c ) e ~ 2 of the integrand in Eq. (|2.10[) at some critical value 
z c . Denoting averages with respect to the new functional 
by (■ ■ -) z , and defining 



g$(z) = ln(exp(gn)) z « q(n) 



(2-11) 



asymptotically, we get by using Jordans lemma that the 
asymptotic behavior, e.g., of V(N) for large ./V is given 
by 



pa+100 . 
V(N) = — exp[zN + g*(z)] 

J a— ioo ^"^ 

(*c) f 

Jo 



z c N+q$(z 



Disc$(z c - x) _ xN 
dx -i— -e 



(2.12) 



where the last row gives the asymptotics for large N. 
Here, z c is the first singularity of <&(z) ~ (z — z c ) e ~ 2 , 
which is a branch point on the negative real axis, and 
the contour of the path integral is deformed into a path 
above and below the branch cut beginning at the singu- 
larity. Disc $ denotes the discontinuity of the function $ 
at the branch cut. The non-universal factor qe z ^ N +i^(. z <^) 
depending exponentially on N is common to all averages 
defined by Eq. (|2.5p and therefore cancels from all mean 
values (O)jv- 

Now, we return to our response functional J z to re- 
fine it into a form that suits us best for our actual field 
theoretic analysis. As discussed above, the gradient term 
proportional to c is redundant. To eliminate this term, 
we apply to the field h the shift and mixing transforma- 
tion 



n(r, t) — > n(r, t) + a — ocrooo(r) , 



(2.13) 



where a is a free parameter at this stage. Ultimately, this 
parameter is defined by (n) = which means that the 
diagrammatic perturbation expansion is free of tadpoles. 
Defining r = r p — ga, p = (g' + fa)a — acr, h = z + r p a — 
ga 2 /2, the stochastic functional J z (|2.9[) takes the form 



Jz 



A A / dth [\- l d t + T-V 2 



gim c 



[a / dth( 
,)n+( 



P 2 , 50 3 
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, 92- 

g-ifn n 



} ■ (2-14) 



Here, we could have set r equal to zero by exploiting 
that a is a free parameter. Instead of doing so, we rather 
keep r in our theory as a small free redundant parameter. 
We will see later on that keeping r comes in handy for 
renormalization purposes. In Eq. (|2.14[) . we have elimi- 
nated couplings that are of more than third order in the 
fields because they are irrelevant. We do not write down 
in detail the relatively uninteresting relations between 
the new third-order coupling constants and the old ones. 
Note that J z contains two similar couplings: g 2 nnm and 
gihnrrioo. Whereas the first coupling respects causal or- 
dering, which means that n is separated by an infinites- 
imal positive time-element from the nw-part resulting 
from the Ito-calculus 25| , the second one respects causal- 
ity only between h and n. In contrast to the m-part, the 
moo-part contains all the n with times that lie in the past 
and in the future of n. This property is the heritage of 
the time-delocalized noise term, and of the introduction 
of the ^-function in Eq. (|2.5[) as a final (t = oo) condition 
that destroyed the causality of J z . Even if we had disre- 
garded the noise term proportional to / in Eq. ()2.2I) ini- 
tially, the fmTOoo-coupling would be generated by coarse 
graining, and hence it must be ultimately incorporated 
into the theory to yield renormalizability. 

The relevance or irrelevance of the different terms in J z 
follows from their dimensions with respect to an inverse 
length scale p, such that time scales as p~ 2 . Fundamen- 
tally, one has to decide which parameters are the crit- 
ical control-parameters going to zero in mean-field the- 



ory. At the collapse transition these are r 



/'" 



and h ~ p^ d+T >l 2 yjj. The dimensions of the fields are 
then given by h ~ m ~ p^ d ~ 2 ^ 2 , and n ~ p( d+2 ^ 2 . It 
follows that all the coupling constants g$, g±, <?2, and 
g' 2 have the same dimension pS % ~ d ^ 2 . Note that n is 
tied always to at least one factor of n as a result of ab- 
sorptivity of the process. Hence, all the terms in J z are 
relevant for d < 6 spatial dimensions, and the upper crit- 
ical dimension of the collapse transition is c? c = 6. The 
situation is different if p is a finite positive quantity, that 
is in the swollen phase. Then p can be absorbed into the 
fields by a rescaling transformation which amounts to for- 
mally setting p = 2. The field dimensions then become 
to - p d ' 2 , n - ^( d + 4 )/ 2 , and h - ^ d ~ 4 )/ 2 . It follows 
that h - p d/2 , go - p~ d/2 , si - ^ (4_d)/2 , 92 - M (2_d)/2 , 
and ^2 ~ /i' 8 ~ d '/ 2 . Hence, in the swollen phase only 
92 = 9 is relevant, now below 8 spatial dimensions. The 
other couplings can be safely removed. Then, J z reduces 
to the dynamical response functional of the usual simple 
GEP. 

Recently, we have shown that g' 2 becomes weakly irrele- 
vant (with a correction exponent of order 0(6 — d)) at the 
RG fixed point corresponding to the collapse transition 
of RBPs Q- Hence, we can neglect this coupling not only 
for the swollen phase but also for the collapse transition, 
even though dimensional analysis suggests its relevance 
for the latter in the vicinity of the Gaussian fixed point. 
In Ref. [5( we argued that hat the vanishing of 172 in - 
dicates that cycles or loops of RBPs are irrelevant. In 
other words: large RBPs in the swollen phase and at the 
collapse transition are dominated by tree-configurations. 

The last step in setting up the response functional for 
the dynamical description of swollen or collapsing RBPs 
is a duality transformation that interchanges absorptivity 
in favor of causality: 

n(r, t) = s(r, — t) , n(r, t) = s(r, —t) , 

/oo 
dt's{r,t'), 

m oc (r) = m(r, 00) =: <^(r). (2-15) 

This step leads to the response functional 

J = f d d xh I dt~s(\- l d t + T-V 2 +g l (p- ! ^s\s 

+ (f^ + f^ + ^)}. ( 2 - 16 ) 

which will serve us in the following as the vantage point 
of our dynamical field theory. 

In this paper, we are mainly interested in the dynam- 
ical aspects of our theory. On occasion, however, we will 
also discuss some of its static aspects. On one hand, 
this will make our presentation more self contained be- 
cause we will need some of the previously derived static 
results as input as we move along. On the other hand, 
we feel that certain elements of the theory are more eas- 
ily discussed statically rather than dynamically. For the 
static aspects, we do not need the full response func- 
tional. Rather it is sufficient to consider the quasi-static 



Hamiltonian 



t> * — =0 



U = 



/ d d xi(p(T - \J 2 )ip + ^if 2 + hip 



90 -3 , ~2 92-2 

—99 +giip ip- —<p>tp 



+ ip(r - V 2 + .91^ - g 2 <p 



)*} 



(2.17) 



that follows from J by setting s(r, t) — > y(r). Here, we 
have added a pair of fcrmionic ghost fields (ipyip) that 
automatically guarantees the original causality rule in 
Feynman diagrams through its couplings to the fields (p 
and tp Q . We note that the quasi-static Hamiltonian has 
BRS supersymmetry which indicates that only tree-like 
polymer configurations are relevant as noted above. 



C. Diagrammatics: SAWs on tree diagrams 

Figure [T] shows the diagrammatic elements of our the- 
ory as resulting from J and %. There is the dynamic 
propagator which reads 



G(q,i)=^)exp[-(T + q 2 )i], 



(2.18) 



in momentum-time representation or, after Fourier trans- 
formation, 



G(q, W ) 



1 



i\uj - 



(2.19) 



in momentum-frequency representation. Its static coun- 
terpart is given by 



G(q)=G(q,0) = 



1 



(2.20) 



The dynamical functional also features the static corre- 
lator 



pC(tf 



~P 



(r + q 2 



(2.21) 



Note the negative sign of the correlator that has its origin 
in the fact that the fluctuations of the fields s and (p are 
purely imaginary (a heritage of the imaginary response 
field h). In addition to these elements represented by 
lines, there are the 3 vertices shown at the bottom of 

Fig.ffl 

Recall that we are primarily interested in scaling ex- 
ponents for transport and diffusion on RBPs. These, we 
will extract from the response function 

Gi,i(r,t) = (n(r,*)n(0,0)} = <s(-r,i)s(0,0)) (2.22) 

and its renormalizations. Its (quasi-)static renormaliza- 
tions are well known from earlier work [4], y, [29|. Its 
dynamical renormalization is known only to 1-loop order 
for the swollen phase as shown by one of us [15J and en- 
tirely unknown for the collapse transition. Hence, it is 



G(q,t) 



q q 

> fw« = 1 — -^--< = C(q) 

q q 



= >-----< =-pC(q) z 
4^P=-lg \ =-Xg, 4—0(=lg 2 

Figure 1. The elements of the dynamic Feynman diagrams 




t<- 



Figure 2. A tree-like subdiagram. 



our main task here to calculate the renormalization factor 
Z\ pertaining to the dynamical coefficient A. To this end, 
we consider the one-particle irreducible amputated self- 
energy diagrams with an outgoing amputated s-leg and 
an ingoing amputated s-leg. Let us forget for the moment 
all the static C-lines of a diagram that are introduced by 
the couplings to the static field (p. Then any of these 
diagrams is reduced to a pure time-ordered tree-diagram 
that has its origin in the outgoing amputated s-leg, see 
Fig. [2j The tree consists of G-lines and ^-vertices. The 
ingoing amputated s-leg is one of the pending s-legs. The 
SAW along G-lines which connects the two amputated 
legs is unique. Note, that the tree-structure with this 
unique SAW is a consequence of the limit g' 2 — > 0, which 
we may hence call the tree-limit. Now we reintroduce 
the G-lines by inserting go-. g\- and p-vertices and satu- 
rating the pending non-amputated s-legs with the static 
<^-legs. Then, integrating over internal times from right 
(earlier times) to left (later times), it is easy to see that 
the internal s-legs are converted by the integrations into 
<p-legs with exception of the s that are part of the con- 
necting SAW between the amputated external two legs. 
The conversion of s-legs into (p-legs turns the correspond- 
ing G-lines into G-lines. In other words, the parts of 
the diagram that do not belong to the connecting SAW 
become purely quasi-static. One can think of the time- 
integrations as having the net effect of decomposing any 
self-energy diagram into a dynamic part, the connecting 
SAW, and a quasi-static residual part. This decomposi- 
tion is visualized in Fig. |3] After the decomposition, the 
frequency dependence of any self-energy diagrams solely 
resides in its connecting SAW. To calculate the dynam- 
ical renormalization factor Z\, we have to calculate the 
parts of the self-energy diagrams that are proportional 
to iuj which we can do, following standard field theoretic 
procedures, by making iwss-inscrtions. After what we 



4 • < — • «— •- 




the response functional reduces to 



or 



Figure 3. Decomposition of a self-energy diagram into a dy- 
namic connecting SAW and a quasi-static part. 



just have learned about the decomposition effect of the 
time integrations, it is clear need to put these insertions 
only into the G-lines forming the connecting SAW. 

All these considerations are easily generalized from 
self-energy diagrams to general one-particle irreducible 
diagram with A* external s-legs and N external s-legs. 
Let us consider such a diagram beginning, as above, by 
temporarily omitting all static C-lines. This leads to a 
decomposition of the diagram into a (spanning) forest of 
A" disconnected time-directed trees rooted in the outgo- 
ing amputated s-legs each featuring a subset of the A" in- 
coming amputated s-legs. Then, Fourier-transformation 
from time-arguments to frequencies shows that the fre- 
quency flowing out through a s-leg is the sum of the 
frequencies flowing into the diagram through the corre- 
sponding set of s-legs. Thus, there are A~ frequency- 
conservation laws. The argument concludes by tying to- 
gether the disconnected trees by readmitting the static 
C-lincs. 



III. THE SWOLLEN RBP 

In this section, we use our model and our insight into 
the structure of the diagrammatic expansion developed 
in Sec. Ill CI to calculate the dynamical exponent z and 
the related transport exponents for the swollen phase. 
Before we embark on this quest, we will briefly review 
some of the quasi-static results that we can utilize as 
input for our dynamical calculation. Furthermore, we 
have performed rational approximations to improve the 
known results for the static exponents that we would like 
to present here along with comparisons to the available 
numerical results for these exponents. 



A. Functionals, renormalizations, and static results 

As discussed above, c/Oj <?ij and g' 2 are strongly ir- 
relevant for the swollen phase, and hence we now set 
go = 9i = g' 2 = 0. Moreover, p is a finite and positive 
quantity that can be reset through a simple re-scaling 
transformation. Hence, we have the freedom to choose 
p = 2 for simplicity and we do so. With these settings, 



U S I 



l d xl 



A / dt§(\- 1 d t + T-V 2 -?-s)s 



{? 



tup)} 



(3.1) 



where g = g 2 .This functional describes besides the dy- 
namical creation of the swollen RBPs the dynamics at the 
Yang-Lee singularity [26[, however with quenched static 
noise \}M L - With the settings for the swollen phase, the 
quasi-static Hamiltonian including ghost-fields becomes 

U sw = / d d x{p>( T -V 2 -^)tp + (p 2 + hip 
+ ${t- V 2 -gtp)i/)\ 
= { d d xd6d6 | -$(r-a)<P--<P 3 + /*$}, (3.2) 

where we have introduced Grassmannian anticommuting 
super-coordinates 9, 9 with integration rules J dd 1 = 0, 
J d91 = 0, J dee = 1, [dee = 1, and defined a super- 
field $(r, 6, 6) = ip(r)+ie-ip(r)+iip(r)6+ee<p(r), as well as 
the super-Laplace operator □ = V 2 + 2dgdg . The Hamil- 
tonian ~H SW shows full supersymmetry, i.e., besides the 
symmetry against super-translations (BRS-symmetry) it 
also has super-rotation symmetry. Parisi and Sourlas 
[4j showed some 30 years ago that the full supersym- 
metry leads to dimensional reduction because it makes 
the Hamiltonian equivalent to the ordinary Yang-Lee- 
Hamiltonian in two lesser dimensions (see also the rigor- 
ous work on the dimensional reduction of Brydges and 
Imbrie [27|, [28[). Exploiting this relation, all static renor- 
malizations are known up to third order |29f . 

We note that in contrast to the quasi-static Hamil- 
tonian which, as we just have seen, can be written in 
the form of the supersymmetric Yang-Lee-Hamiltonian 
by the introduction of a superfield, the dynamic response 
functional is not supersymmetric. Hence, we unfortu- 
nately cannot exploit dimensional reduction in our dy- 
namic calculation. 

In the following, we use dimensional rcgularization. 
For the swollen phase, we employ the renormalization 
scheme 



s^s = Z 1/2 s, S^ 1 = Z X Z 1/2 S, 


(3.3a) 


X^k= Z^X, T^T = Z~ 1 Z t t, 


(3.3b) 


g^g = Z- 3 / 2 Z g g, G e g 2 = utf , 


(3.3c) 


h^h = Z- 1 ' 2 (h- 9 ~Ai 3 ), 


(3.3d) 


(p>, ip, $, ip) -> (<p, p>, $, -0) = z 1/2 (<p>, ip, $, vo . 


(3.3e) 



Not all the renormalization Z-factors in the renormal- 
ization scheme are independent. The form-invariance of 
J sw and Ji sw under the shift s — > s + a leads to the Ward 
identities Z g = Z T and Z T = 1 + g 2 A [2g|. As mentioned 
above, the Z-factors other than Z\ are known to 3-loop 



order. In our dynamical calculation 2-loop calculation, 
we will exploit the known result for the field renormal- 
ization Z to 2-loop odcr as an input 



Z 



u 
3e 



^)©W>, 



(3.4) 



where e — 8 — d. 

For completeness, let us mention here that the known 
3-loop result for the entropic scaling exponent 9 as fea- 
tured in Eq. (|2~%|) reads $2% 
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(3.5) 

The scaling exponent Up of the gyration radius of the 
branched polymer, Rn ~ N VF , is related to 9 by 







9=\ 


[d - 2)up 


+ 1, 




(3.6) 


The 


resulting e 


-expansion of vp : 


reads 






Up = 


1 e 


29 2 


+ /C0) 

i486 


8879 \ 3 
7558272/ 


+ 0(£ 4 ). 

(3.7) 


4 + 36 + 


23328° 



Now, we improve these results by making rational ap- 
proximations that incorporate the exact results that are 
known for 9 in dimensions d from 1 to 4, namely 9{d = 
1) = 0, 9{d = 2) = 1, 9{d = 3) = 3/2, 0(d = 4) '= 11/6. 
We get 

5 e / 1 + 1.37622e - 0.0130833e 2 - 0.0171653e 3 \ 

^ 2 ^ 12 V 1 + 1.13239e - 0.210607e 2 + 0.00685362e 3 / ' 

(3.8) 
For up, we include the exact values up{d = 1) = 1, 
u P (d = 3) = 1/2, up(d = 4) = 5/12 and get 

_1 e_/ l + 1.17714e-0.113178e 2 \ 

^ P ~ 4 + 36 V 1 + 1.13239e - 0.210607e 2 + 0.00685362e 3 / 

(3.9) 
With regard to the relation (|3.6p , we note that our inde- 
pendent rational approximations satisfy these relations 
to order O(10 -6 ). In Table HI we compare the numerical 
values resulting from these approximations for various di- 
mensions to recent simulation results by Hsu, Nadler, and 
Grassberger [9|, [jjj] . Over all, the two agree remarkably 
well over a wide range of dimensions. 



B. Calculation of the dynamical scaling exponent 

Now we return to our main goal, the 2-loop calculation 
of the dynamical exponent z via the calculation of Z\. 
As a warmup, let us first do a quick 1-loop calculation 
that reproduces the known 1-loop result. At this order, 
there is only one self-energy diagram, namely diagram 
(1) of Fig. (|4]) . Using the graphical elements pictured in 
Fig. (JTJ), this diagram translates into the formula 



(1) - Ag 2 



(-2) 



p (r + p 2 ) (iu/X + t + (p - q) 2 



(3.10) 




Figure 4. Selfenergy diagrams for swollen RBPs to 2-loop 
order. 



where f ... is an abbreviation for the d-dimcnsional mo- 
J p 

mentum integral (27r)~ d J d d p .... The calculation of this 
and in particular that of the higher order diagrams can be 
simplified be resorting to massless propagators. Practi- 
cally, this is done by expanding the integrands in powers 
oiiut/X and r. For the 1-loop diagram (1) of Fig. (0} this 
gives 



(1) = -2A 5 2 



iu/X + t 



'p l(p-q) 2 (p 



2 n 2) 



((p-q) 2 ) (P 2 )' 



2r 



(p-q) 2 (p 2 ) 



(3.11) 



Note that the massless integrals produced by the expan- 
sion are not infrared divergent as long as d > 6. For the 
calculation of these integrals, we refer to Appendix [A] 
Their e-expansions about d = 8 lead to the renormalized 
vertex function to 1-loop order 

r^ L) (g, u) = (ZZxilu + Z t Xt + ZXq 2 ) 

- ^- (2iu + 6Ar + Xq 2 ) + . . . . (3.12) 
3e 

Hence, we obtain to this order 

Z = Z X = 1 + ^-, z T = l + — , (3.13) 

3e e 

which leads to the 1-loop result z — 2 + r\ + 0(e 2 ) (15I ). 
where r\ = — e/9 + 0(e 2 ) is the anomalous dimension of 
the field s HI. 

Now turn to the 2-loop part of the calculation for which 
we consider the 2-loop diagrams of Fig. (TTJ). The parts of 
these diagrams linear in iui are 

(2a) + (26) + (2c)) = (2F(5, 2, 1) + ^F(4, 2, 2) 

x (-Aiujg 4 ) , (3.14a) 
(2d)u, = (2F(5, 2, 1) + F(4, 2, 2)) (-Aicog 4 ) , (3.14b) 
(2e) lu = (2F(4, 3, 1) + F(3, 3, 2)) (-Aicog 4 ) , (3.14c) 
((2/) + (2.g) + (2/0) ^ - (2F(4, 3, 1) + 2F(4, 2, 2) 

+ 2F('3,3,2))(-4iLog 4 ), (3.14d) 



Table I. Rational approximation estimates of the exponents 9 und vp compared with simulation results by Hsu, Nadler, and 
Grassberger [jj. 



d 


1 


2 


3 


4 


5 


6 


7 


8 


e 





1 


3/2 


11/6 


2.0769 


2.2603 


2.3986 


5/2 


Hsu et al. 








1.835(6) 


2.080(7) 


2.261(12) 


2.40(2) 




Vp 


1 


0.64267 


1/2 


5/12 


0.35896 


0.31507 


0.27973 


1/4 


Hsu et al. 




0.6412(5) 




0.4163(30) 


0.359(4) 


0.315(4) 


0.282(5) 





where the F(. . .) are the frames of the diagrams, i.e., 
their parts consisting only of the momentum integrations 
without any coupling constants or symmetry factors. For 
a more precise definition of the F(. . .), we refer to Ap- 
pendix [AJ Next, we calculate the counter-terms of these 
frames. For the simplicity of the argument, let us just 
say here that we apply some calculation procedure £ to 
the frames that produces as its result the counter-terms 
of the frames: 



CF(...) = C(. 



(3.15) 



A precise definition of €. along with some details of the 
calculations is given in Appendix [3] Application of £ 
results in 

7 



4(7(5, 2,1) + -(7(4, 2, 2) 



c((2o) + • • • + (2fc) 

+ 4(7(4, 3, 1) + 3(7(3, 3, 2)) (-4iw<? 4 ) 
43 



9e 2 V 132* 



e )iuju 



(3.16) 



These counter-terms yield the 2-loop contribution to the 
renormalization-factor product 



ZZ 



2u 11 



43 



- 1+ S + ^( 1 -^)" 3+0 <" 3 >- <3 - 17 » 



Using Eq. (|3.4[) . we finally obtain the wanted dynamic 
rcnormalization factor 

z ^ 1 + 5 + i( 1 -f)" 2+0< " 3 ^ <3 - 18) 

Next, we discuss the renormalization group equation 
for the vertex functions and its solution. As the result of 
the dynamic tree structure of the Feynman diagrams, the 
vertex functions r^ N ({u>/\}, {q},r, u, fj) depend only 
on the incoming frequencies of the N s-legs. Setting these 
frequencies to zero, we get in the case N = 1 the static 
vertex functions 



r^; t) ({q}) = r ljW ({ w /A = o},{q}) 



(3.19) 



that are related by dimensional reduction to the vertex 
functions r) v+1 ' of the Yang-Lee theory in two lesser di- 
mensions. The RGE of the dynamical T^ N reads [26| 



U V L N,N — I o N,N 



7r^-<5iV,l^JV,0, (3.20) 



where 

£V = A*fy + CA5 A + f3d u + KTd T (3.21) 

is the RG differential operator, and 

7 = ndp In Z\ , 7 T = /j<9 m In Z T \ Q , C = M, In Z x \ , 
K = 7-7r, /3 = (-e + 37-2 7T )u. (3.22) 

For the Wilson functions featured here, we have the ex- 
pansions 

7 = -| + ^u 2 + 0(u 3 ) , 7t = -2u + — w 2 + 0(u 3 ) , 



C = -| + ^ 2 + 0( U 3 ) 



(3.23) 



The fixed point u* of the RGE determined by /3(u«) = 
reads 



= * + tS £2 + °( £3 ) • 



(3.24) 
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At this fixed point, it follows that 

e — 2k(w„) = 7(w*) = 77. (3.25) 

Hence, the static exponent 77 and the dynamic exponent 

z = 2 + C(u*) (3.26) 

are the only independent critical exponents of the prob- 
lem. 

Shift invariant observables free of redundancies are 

M = T-g{s), H = 2gh + r 2 , (3.27) 

with RGEs 

T> i1 M = kM, D M ff = (7 + k)# • (3.28) 

Note that M and H are linearly related to the Laplace 
transform $(z) of the cluster probability V(N) and the 
original Laplace variable z, respectively. M replaces the 
redundant parameter r in all shift invariant quantities. 
The integration of the RGEs (|3.28[) at the fixed point 
yields the equation of state [26[ 



M = const. H 



0-2 



(3.29) 



Likewise, by integrating the RGE p.20j) at the fixed 
point, we obtain 

F(tlr z rH' /p ) 
(«(r,t)5(O,0))= { ' ' j (3.30) 



r 



d-2+ri+z 



Table II. Rational-approximation estimates of the fractal dimensions of swollen RBPs compared with numerical results of 
Havlin, Djordjevic, Majid, Stanley, and Weiss [g). 



d 


1 


2 


3 


4 


5 


6 


7 


8 


(l w 


2 


2.7138 


3.3127 


3.8638 


4.3960 


4.9238 


5.4560 


6 


Havlin et al. 




2.78(8) 


3.37(10) 


3.89(12) 










flmin 


1 


1.1578 


1.3127 


1.4638 


1.6101 


1.7499 


1.8811 


2 


Havlin et al. 




1.17(8) 


1.36(10) 


1.49(12) 










df 


1 


1.5560 


2 


12/5 


2.7859 


3.1739 


3.5749 


4 



■D 1. 



1.6 



1.4 



o 1.2 



12 3 4 5 6 7 8 
spatial dimension d 

Figure 5. (Color online) Rational-approximation estimate of 
the fractal dimension of the minimal path compared with nu- 
merical results of Ref. [g|. 




3 4 5 6 7 
spatial dimension d 

Figure 6. (Color online) Rational-approximation estimate of 
the fractal dimension of the random walk compared with nu- 
merical results of Ref. [g. 



for the response function. The static exponents are re- 
lated by 



Vp 



-2 = 



'/ 



d-A + ri 
d — i] 



(3.31) 



in conformity with the relation (|3.6|) . For the dynamical 
exponent, we obtain the 2-loop result 



z = 2 



«->= 2 -f-§(!) 2 +°<^ t" 2 ' 



Finally, we improve the numerical accuracy of this re- 
sult by incorporating the exact value z = 1 for d = 1 
(e = 7) through the rational approximation 



1134 + 245e 
10206 +1391e 



(3.33) 



Using relation (|1.4|) together with df = 1/fp and in 
conjunction with the rational approximations given in 
Eqs. (|3.9[) and ()3.33[) . we obtain the fractal dimensions 
displayed in Table UU and Figures [5] and El These figures 
also show the available simulation results [6J. We find 
that our field theoretical results agree remarkably well 
with the latter. 



IV. THE COLLAPSING RBP 

Now, we turn to the collapse transition. As we did for 
the swollen phase, we will first briefly review some known 



results |5| that we need as input as we move along. The 
main task will be once again the calculation of the dy- 
namical exponent z via the calculation of the dynamical 
renormalization Z\. 



A. Static renormalizations and results 

As discussed in Sec. UU g$, g\ and c/2 are all relevant 
for the collapse transition and hence need to be kept in 
the field theoretic functionals. As also discussed, g' 2 is 
weakly irrelevant at the collapse transition, and hence 
we can take the tree limit g' 2 —> here. Thus, we work 
with the response functional J as given in Eq. (|2.16[) 
and, as far as the static properties are concerned, the 
quasi-static Hamiltonian Ti. as given in Eq. (J2.17I) . 

Our renormalization scheme for the latter in dimen- 
sional rcgularization is 

{ip, ip, xjj, i/0 ->■ (tp, ip, i/5, i>) = Z 1/2 (if, ip + Kip, -0, ip) , 

(4.1a) 
Z^L= Z^g-L, (4.1b) 

h -> h = Z- l ' 2 (h + ici/ V £/2 Z • A ■ L) , 

(4.1c) 
G l J 2 g a -> G\/ 2 g a = Z-*' 2 {u a + B Q ) M e/2 , (4.1d) 

where r = (p, r) and where e now measures the deviation 
from d = 6, e = 6 — d . Since we are not interested 
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here in the renormalization of the control parameter p 
that defines the cross-over variable to the swollen phase 
and which goes to zero at the collapse- or 6>-linc, we set 
p = in the following calculations. Scaling invariant 
combinations of the coupling constants are defined by 



V = UiU 2 , 



u ul, 



(4.2) 



with fixed point values 

u, = 0.6567 (e/6) + 2.9707 (e/6) 2 + 0(e 3 ) , (4.3a) 



u;, = 0.7052 (e/6) 2 + C>(e 3 ) 



(4.3b) 



As it did for the swollen phase, the shift-invariance leads 
to Ward identities for the collapse transition. Here, how- 
ever, these Ward identities do not result in a scaling re- 
lation between the polymer exponents 9 and vp but they 
nevertheless simplify the calculations or provide consis- 
tency checks. The e-expansions for these exponents are 



9= - - 0.4925 (e/6) - 0.5778 (e/6) 2 , (4.4a) 



v P = - + 0.1915 (e/6) + 0.0841 (e/6)" 



(4.4b) 



B. Calculation of the dynamical scaling exponent 
of collapsing RBPs 

For our dynamical calculation, we complete the renor- 
malization scheme (|4.1|) by setting 

s^s = Z 1/2 (s + K<p), (4.5a) 

s -> s = Z x Z 1/2 s , X^X = Z^X. (4.5b) 

Let us first describe the 1-loop part of our calculation. 
The diagram (1) of Fig. (|7J leads to the following contri- 
bution to the sclfencrgy expanded to first order in u and 

T 

1 




p (^/A + T + p2)(T+(q- P )2) 

1 ico/X + t 

P i (p - q)V ( P -q) 2 (p 2 ) 2 



p 2 ((p-q) 



(4.6) 



The factor 2 stems from the two possible orientations of 
the propagator in diagram (1). Performing the integra- 
tions as described in Appendix [B] and carrying out an 
e-expansion about d = 6 we find the renormalized 1-loop 
vertex function 



T\y(q, u) = (ZZ x iuj + Z t Xt + ZXq 2 

- —(3iuj + 6XT + Xq 2 ) + 
3e v 



(4.7) 




Figure 7. Diagrams to 2-loop order for the collapse transition. 
Dotted lines symbolize all possible arrangements of the static 
propagator including ghosts. 



Form this, we read off the 1-loop renormalization factors 

2v Av Av 

Z = l + — , Z T = 1 + — , Z x = l + — , 4.8) 

3e e 3e 

with Z and Z T as already known from [5(. 

Now we turn to the 2-loop self-energy diagrams of 
Fig. [3 For details of the calculation, we refer to Ap- 
pendix [B] It encompasses the frame counter-terms shown 
in Fig. (|12D for which we get the results compiled in 
Eq. (|B3|) . From these, we obtain the following results 
for the counter-terms of the Fcynman diagrams: 

C(2a) c = 4A( 5l52 ) 2 {Ci(3,l,l) 

- (2Ct(3,2,l) + 3C*i(4,l,l))T 

- (Ct(3, 2, 1) + 2Ct(4, 1, l))y } , (4.9a) 
C(2b) c = X(7( gi g 2 ) 2 -. 90 .92 2 ){Ci(3,l,l) 

- (2Ct(3, 2, 1) + 3Ct(4, 1, l))r - Ct(4, 1, l)y } , 

(4.9b) 

C(2c) c = A(9(. 91 . 92 ) 2 - . 90 .92 2 ){ci(2,2, 1) 

- (4Ct(3, 2, 1) + Ct(2, 2, 2))r - 2Ct(3, 2, 1)^} , 

(4.9c) 

C(2d) c = AX{g l92 ) 2 \ct{2,2,l) 

- (4Ct(3,2,l) + Ct(2,2,2))r 

- (2Ct(3, 2, 1) + Ct(2, 2, 2)) y } . (4.9d) 



These add up 


to 








C((2a) + 


■■■ + (2d)) c 






= 


v67 
A 18 


191 \ 2 

e)v - 

216 J 


■(1 


13 \ i 

e w 

216 / J 




y63 

At - 


301 \ 2 
- -77r£)v - 
24 / 


(!- 


23 \ i 

tttE w 

24 / . 


Xt 


+ 


T 2 " 




-/25 

At" 


103 \ 2 
24 / 


(I- 


11 X " 

— e ]w 
72 J - 


iuj 


+ 


e^ 



M 2 



(4.10) 



This defines the 2-loop contributions to the renormaliza- 
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Table III. Estimates of the fractal dimensions of collapsing 
RBPs obtained via rational approximations. 



d 


1 


2 


3 


4 


5 


6 


CL W 


2 


3.038 


3.882 


4.631 


5.332 


6 


(*mi n 


1 


1.192 


1.396 


1.612 


1.824 


2 


df 


1 


1.846 


2.486 


3.019 


3.508 


4 



tion factors Z, Z T , and ZZ\. Finally we obtain 



Z=l + — + ^r 



2v 1 
3e e 2 


/67 
A18~ 


191 \ 2 
216 / 


/ 5 13 ^ 

~ V18 ~ 216 £ / 

(4 


.11 


4v 1 

£ £ 


/63 
AT" 


301 \ 9 
- -— -e hr 
24 / 


/5 23 \ 

" 12 ^ 24 e r 

(4 


) 

.11 


4v 1 
3e e 2 


v83 
\~9~ 


92 \ 2 
27 / 


(|-^ e ) W . 





^r = l + — + ^ hr-177 



(4.11c) 

As they should be, the first two are in conformity with re- 
sults of [5| . The dynamic renormalization group function 
becomes 



c = 



dkiZx 



d In// 



4v 

y 



/184 2 5 
V 27 12 



0(3-loop) . 

(4.12) 

Using the fixed point values of the coupling constant 
given in Eq. ([4.3(1 , we obtain the e-expansion of the dy- 
namic exponent and by the same token the scaling di- 
mension of the minimal path: 

rfmm = z= 2+C* = 2-0.8756 (e/6)-1.3162 (e/6) 2 +0(e 3 ) 

(4.13) 
Taking into account the exact values z = vp = 1 for 
d = 1 (e = 5), we propose the rational approximations 



^min — 



2- 0.14593 e 



l + 0.726717e 
1 + 0.476179 e' 



and 



\/df = Vp « - + - 



0.0319167 e 



0.0733194e- 0.016825 e 2 



(4.14) 



(4.15) 



the latter being based on the e-expansion of Up as given 
in Eq. (|4.4[l . Table IIIII lists the numerical values result- 
ing from these rational approximations for various spatial 
dimension. 



V. CONCLUDING REMARKS 

In summary, we have studied diffusion and trans- 
port on swollen and collapsing randomly polymers us- 
ing renormalizcd dynamical field theory. In particular, 
we have calculated the diffusion exponent and a set of 
transport exponents including the fractal dimension of 
the minimal path and the resistance exponent to 2-loop 



order. For the swollen polymer, our results are an im- 
provement of the previously known 1-loop results. For 
the collapse transition, our results are entirely new in the 
sense that hitherto no results beyond mean-field theory 
existed. 

From a conceptual or diagrammatic standpoint, it was 
interesting to sec that the dynamical Fcynman diagrams 
for the self-energy decompose into dynamic part which 
has the form of a SAW connecting the external legs and 
a residual quasi-static part. This observation simplified 
the calculation of the dynamical renormalization enough 
to enable us to treat the problem to 2-loop order. 

As far as we know, the existing numerical results for 
diffusion and transport on substrates in the universality 
class randomly branched polymers and lattice animals 
are more than 25 years old. For the collapse transition, 
no such results exist at all. We hope that our results 
encourage simulation work leading to improved and ex- 
tended numerical results. 
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Appendix A: 2-Loop massless calculation of frames 
and counter-terms near d — 8 

Here, we calculate the frames of the dynamical dia- 
grams for the swollen phase using dimensional regular- 
ization [30(. We consider only those frames that are re- 
quired to determine the dynamic renormalization factor 
to 2-loop order. The calculation can be greatly simplified 
via the consequent use of massless propagators [31( . Of 
course, massless propagators bring about the danger of 
infrared (IR) singularities which would require the intro- 
duction of IR-counter-terms [32|, |33( if they indeed oc- 
curred. However, cleverly using infrared rearrangements 
of external momenta [34|, [35| , the IR-singularities can be 
avoided at least up to 3-loop order. Figure [8] lists the 
massless frames that we need for our 2-loop calculation. 
The wiggly lines symbolize external momenta ±q flowing 
into or out of the frame. These wiggly lines are positioned 
so that the frames are free of IR-singularities. 

First, let us return to the 1-loop part of our calculation 
sketched in Sec. IIIII Equation p. lip for diagram (1) can 
be expressed as 

(l) = 2A.g 2 {-G d (l,2) + G d (2,2)(^/A + r) 

+ G? d (l,3)2r +...}, (Al) 

where Gd(- • • ) is defined through the fundamental inte- 
gral 

, 2^,/ ^a, - <fc(Ax, W-^+*> . (A2) 

p (p 2 ) ((p-q) 2 ) 



12 



F(3) 



F(4) 



F'(4) 





t' 



F(5,2,l) F(4,3,l) F(4,2,2) F(3,3,2) 



Figure 8. Frames of the diagrams for swollen RBPs up to 
2-loop order. 



This integral can be easily carried out using standard 
methods with the result 



Gd(Ai, A2) 



r(d/2 - Ai)r(d/2 - A 2 )r(Ai + a 2 - d/2) 



(4^/ 2 r(A 1 )r(A 2 )r(d-A 1 -A 2 ) 



e-expansion in 



d leads to 



(A3) 



F(3) = G - (l,2)g a - = -5 a ^ 



( 1+ r 



G e ^ 7 
(A4a) 



F(4) = G„(l,%- e = Ml!(l + g £ + ...)) G eM - 



(A4b) 



F'(4) = G d (2, 2)g" £ = ^1 (l + jj e + . . .)) G e ^ , 

(A4c) 

with G e = T(l + e/2)/(47r) d / 2 . Note that for these values 
of the parameters Xi , the integral (|A2[) does not contain 
IR-singularities. Hence, the e-poles in Eqs. (|A4j) purely 
arise from UV-singularities as they should, i.e., we have 
successfully avoided IR-singularities even though we re- 
sorted to massless propagators. 

Next, we consider the 2-loop diagrams. The frames 
with the arrangement of external momenta as shown in 
Fig. © are free of IR-singularities. They are easily inte- 
grated through a successive application of Eq. (|A2[) . We 
obtain 



F(5, 2, 1) = G d (2, 6 - d/2)G d (l, 2)<T 2£ 

)g 2 m _2£ , (A5a) 



(/V<?) 2e / 25 
36e 2 V 12 



F(4, 3, 1) = G rf (2, 6 - d/2)G d (l, 3)<T 2£ 

_ ( M /g) 2g / 11 n 2 2e 

F(4, 2, 2) = G d (2, 6 - d/2)G d (2, 2)q 



(A5b) 



-2_^ 



0*/«) 



2? 
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lg£2 v l + - £ + ...jG^-, (A5c) 

2e 



F(3,3,2) = G d (l,7-d/2)G d (2,3) g 

_ ( M /g) 2e 



24e 



1 



2, ,-2e 



Gin 



(A5d) 



C(3) C(4) 

Figure 9. 1-loop counter-terms in d — 8 — e dimensions. 



Now, we calculate the counter-terms using a BPHZ- 
type of renormalization [3J, l36| . i.e., a recursive renor- 
malization procedure defined by a sequence of operative 
steps. To this end, we define a pole-separating procedure 
*P operating on e-expansions 



*E 



CiS = 



E 



CiS 



(A6) 



The counter-terms are constructed by the operation 

€ = *&# , (A7) 

where the incomplete renormalization 9V operates on the 
momentum- integrals It of the one-particle irreducible di- 
agrams r. It is recursively defined by 



tfl r =lr+ J2 J r/{7}- II ("^ 

{y}eD r 7.e{7> 



(A8) 



Here Dp is the set of all collections {7} = {71,72, • ■ •} 
of disjunct superficially divergent one-particle irreducible 
subdiagrams 7$ of T, and r/{7J is the diagram obtained 
from r by collapsing each of these subdiagrams to points. 
Applying d\' to a superficially convergent diagram T, 
Dt'ir is finite for e — > 0, and £ir = 0. The operation 
£R' applied to superficially divergent diagrams produces 
poles in e with coefficients polynomial in the external 
momenta and internal masses without non-primitive log- 
arithmic terms. Finally, the operation C/r yields the 
wanted counter-terms, and the complete renormalization 
DUr = SH'ir - £ir = (1 - 93)9*' ir leads to a finite result. 
Of course, applying 9t' to 1-loop diagrams is trivial. 
The counter-terms (dashed boxes) are shown in Fig. [9] 
Using Eqs. (|A4[) . we get 



G(3) = <£F(3) 



G e /z- £ 2 

"9 : 



(if 



G(4) = qjF(4) = «PF'(4) 



G E //- 
3e 



(A9a) 
(A9b) 



Note that we leave factors G e yT e un-cxpanded because 
they arc absorbed by the dimension-bearing coupling 
constants, cf. the renormalization scheme (|3.3[) . Note 
also that because of the application of the operation €, 
it does not matter at which vertices we have injected 
the external momenta as IR-regulators into a given su- 
perficially logarithmically divergent diagram T as long as 
IR-divcrgencies are avoided. 

At 2-loop order, the procedure *p produces the counter- 
terms listed in Fig. [TUJ Mathematically, these stand for 
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C(5,2,l) 



C(4,3,l) 



C(4,2,2) C(3,3,2) 

Figure 10. 2-loop counter-terms in d = 8 — e dimensions. 



C(5 J 2,l)=q}[F(5,2,l)-G d (2,2)C(3)] 

m r (p/q) 2e f, , 25 \ (ji/gy / 5 A\^ -^v 



(Gef*- 



36e 2 



5 

1 -T2 £ 



(AlOa) 



C(4, 3, 1) = q3[F(4, 3, 1) - G d (2, 2)C(4)] 



V{ 



18e 2 I 6 ' 



(G e/ u- 



18e 2 



'-!« 



k> ( 1+ r)}< G -"-) 2 

(AlOb) 



C(4, 2, 2) = qj[F(4, 2, 2) - G d (2, 2)C(4)] 



M A 1 



18e 2 



12 



C(3,3,2)=q3[F(3,3,2)] = 



(G eM - 



24e 



(AlOc) 



(AlOd) 



Appendix B: 2-Loop calculation and counter-terms 
near d — 6 



The counter-terms needed for the collapse transition, 
see Figs. [TT] and [12l can be obtained through a masslcss 
calculation using methods similar to those explained in 
Appendix |A1 in co njun ction with t'Hooft and Veltman's 
'partial p' method [30(. Partially, they can also be ex- 
tracted from results of de Alcantara Bonfim et al. [29[ 
for the usual </> 3 -ficld theory. 

Let us start here by revisiting Eq. (|4.6[) . In terms of 
the fundamental integral (|A2[) . this 1-loop diagram can 



Ct(2) Ct(3) 

Figure 11. 1-loop counter-terms in d = 6 — e dimensions. 




Ct(2,2,l) 



Ct(2,2,2) 




Ct(3,l,l) 



Ct(4,l,l) 



Ct(3,2,l) 
Figure 12. 2-loop counter-terms in d — 6 — e dimensions. 

be expressed as 

(l) c = 2X 9l g 2 { - G d (l, 1) + G d (2, l)(*w/A + 2r) + . . . } . 

(Bl) 

After e-expansion about d = 6, we extract the counter- 
terms 



Ct(2) 



G e fJ, e 2 /~ij.fr>\ G e [i 



-q, Ct(3) = 



3e e 

For the 2-loop counter-terms, we obtain 

.2 



Ci(2,2,l) 



[Ge^ 



3e 2 



Ct(3, 1,1) = - 



(G efi - 



18e 2 



C*i(4,l,l) = 



Ci(3,2,l) 



Ct(2,2,2) = 



{Gefi- 



6e 2 
(G^~ E 
2e 2 
{G sl x-i' 



12 r 



(B2) 

(B3a) 
(B3b) 
(B3c) 
(B3d) 
(B3e) 



Appendix C: Relations between Feynman integrals 

in d and (d — 2) dimensions 

It is well known that supersymmetry relates many 
Feynman integrals with insertions to the corresponding 
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integrals in 2 lesser dimensions [J,[37|- Therefore, it arises 
the question whether such relations can be used to de- 
termine the counter-terms of the frames needed in our 
dynamical calculation. We will show explicitly that this 
indeed is the case at least to 2-loop order. This observa- 
tion may open the route to a 3-loop calculation since all 
counterterms near d = 6 are known [29| . For background 
information on the following reasoning, see the textbook 
by Itzykson and Zuber |36| . 

Consider a 1-particle irreducible Feynman diagram G 
without tadpoles consisting of V vertices, / internal lines, 
and L = I — V+l loops in d dimensions. Each line I car- 
ries a propagator l/(r; + q/ 2 ) where q/ is a d-dimensional 
momentum vector, and 77 an auxiliary mass squared 
which is finally set to zero in a massless calculation. Next, 
let us introduce an (arbitrary) orientation of each line, 
and we define the incidence matrix (e v i), 

1 if the vertex v is the starting point of line I 
£ v i = \ — 1 if the vertex v is the endpoint of line I 

if I is not incident on v . 

(CI) 
Momentum conservation at each vertex v is expressed as 



^Evim = Qi 



(C2) 



1=1 



where Q„ is the external momentum flowing into ver- 
tex v. Only (V — 1) of the conservation laws (|C2[) are 
independent since 



E e ^ = °' (C3) 

The last one produces the overall conservation 

v 
^Q„ = 0. (C4) 



The contribution I(G, {Q„})d of the diagram G to 
the corresponding vertex function, can be written in a 
Schwinger parametric form as [36[ 



si 



I(G, {Q v }) d = rt[d Sl [ cxp ( - V(r, + <£ 

JO l=1 J{q} V , 

V-l 

x ][ ((2n) d 5W(£ £ «*il-Q*))> ( C5 ) 



v=l 



where we focus on the diagram's frame (integral over the 
loop-momenta) and omit any symmetry factors and cou- 
pling constants. Using the integral representation of the 
(^-functions and performing all the arising Gaussian inte- 



gration we arrive at 



I(G,{Q v }) d = / l(c- T ' s 'ds 
Jo 1=1 



exp -Q({s,},{Q„})/P({ a ,}) 



((47r)^({ S/ })) d/2 
Here, V({si}) is given by 

V({s 1 }) = (Att)- l Y[(4tts 1 ) 



(C6) 



1=1 



(4tt) 



v-i 



(C7) 



A is the (V — 1) X (V — 1) matrix with elements A vv / = 
J2i Svis^Sy'i where v,v' = 1, . . . , V - 1. Q{{si}, {Q v }) 
is the bilinear form constructed from the {Qi,} and the 
inverse of the matrix A: 



Q({sih{Qv}) = v({si}) 



V-l 

E 



Qu • (A )„i/Qi/ • (C8) 



V and Q are readily simplified using Kirchhoff 's laws on 
linear electrical networks 1381. 



p({«i}) = ETT 



si 



g({«imq„})= e n 



T l£T 

«(E 

(Ti,T 2 )U0Ti,Ta ^GTi 



Q,)' 



(C9a) 

(C9b) 



Here, the sum T runs over all spanning trees, and the 
sum (Ti, T2) runs over all pairs of mutually disconnecteds 
panning trees (a forest with two trees) on the diagram G, 
respectively. Note that the dependence of I(G, {Q v })d as 
given in Eq. (|C5[) on the dimensionality d entirely rests 
in the exponent of the denominator. 
Now, let us define the procedure 



r = (4ir) L T({-d/dn}) 



(CIO) 



that produces the sum over all diagrams obtained by in- 
serting s 2 into all lines of G not belonging to any tree. 
Applying this procedure to I(G, {Qu})d transforms this 
integral into its (d — 2)-dimcnsional counterpart: 



T*I(G,{Q v }) d = I(T*G,{Qv})d 
= /(G,{QJ) d _ 2 . 



(Cll) 



Next, we turn to the self-energy diagrams in d dimen- 
sions. These diagrams have only 2 external momenta, 
{Qv} = {q, — q}, and the bilinear form Q({si},{Q v }) 
can contain only those partitions into two disconnected 
trees (T\,T2) that disconnect the two external vertices 
where the external momenta enter or exit the diagram. 
The product si runs therefore over the L + 1 lines 

that belong to a cut-set which divides the diagram into 
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two tree-like subdiagrams each containing one external 
vertex. In this case, the form Q({si}, {Q„}) reduces to 

Q({si},{Qv})= £ J] W= : 2*(M)q 2 : 

(C12) 
where v\ £ T\ and V2 G T2. Now, we define the procedure 



C* = {^) L Q*{{-d/d Tl }) 



(C13) 



which produces the sum over all diagrams obtained from 
G by placing an insertion into any of its cut-set lines. 
Applying this procedure to I(G, {q, — q})d results in 

C*/(G,{q,-q}) d = /(C*G,{q,-q}) d 
d 



dq 2 



7(G,{q,-q}) d _ 2 . (C14) 



Note that this theorem relates quadratically diverging di- 
agrams to logarithmic diverging diagrams in 2 dimensions 
higher. 

Having derived the theorems (|C11J) and (|C14|) , we now 
can use them to extract relations between the counter- 



terms encountered in our 2-loop calculations. Applying 
the T* -procedure, we obtain 



Ct(2,2,2) = 12G(3,3,2)), 

Gt(3, 2, 1) = 6G(4, 3, 1) + 3G(4, 2, 2) + 2G(3, 3, 2) , 

Gt(4, 1, 1) = 8G(5, 2, 1) + G(4, 2, 2 . (C15) 

The C* -procedure provides us with the relations 



_d 
dq 
_d_ 

()q 



2 Ct(2, 2, 1) = 4G(4, 3, 1) + 2G(4, 2, 2) + 2G(3, 3, 2) , 



2 Ct(3,l,l)=4G(5,2,l) + G(4,2,2) 



(C16) 



These five relations determine all the four 2-loop counter- 
terms near 8 dimensions from the five 2-loop counter- 
terms near 6 dimensions and in addition to that provide 
us one supplementary consistency check for our calcula- 
tions. 
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